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Abstract 

The mycotoxin aflatoxin is a secondary metabolite and potent human carcino- 
gen. We investigated one mechanism that links stress response with coordinate 
activation of genes involved in aflatoxin biosynthesis in Aspergillus parasiticus. 
Electrophoretic mobility shift assays demonstrated that AtfB, a basic leucine 
zipper (bZIP) transcription factor, is a master co-regulator that binds promot- 
ers of early (fas-l), middle (ver-1), and late (omtA) aflatoxin biosynthetic genes 
as well as stress-response genes (mycelia-specific catl and mitochondria-specific 
Mn sod) at cAMP response element motifs. A novel conserved motif 5'-T/GNT/ 
CAAG CCNNG/AA/GC/ANT/C-3' was identified in promoters of the aflatoxin 
biosynthetic and stress-response genes. A search for transcription factors identi- 
fied SrrA as a transcription factor that could bind to the motif Moreover, we 
also identified a STRE motif (5'-CCCCT-3') in promoters of aflatoxin biosyn- 
thetic and stress-response genes, and competition EMSA suggested that MsnA 
binds to this motif Our study for the first time provides strong evidence to 
suggest that at least four transcription factors (AtfB, SrrA, AP-1, and MsnA) 
participate in a regulatory network that induces aflatoxin biosynthesis as part of 
the cellular response to oxidative stress in A. parasiticus. 
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Introduction 

Reactive oxygen species such as singlet oxygen ('O2), super- 
oxide (©2^), hydrogen peroxide (H2O2), and hydroxyl rad- 
ical (HO') cause oxidative damage to protein, DNA, and 
lipids through protein inactivation and lipid peroxidation 
in cells. Reactive oxygen species (ROS) is produced either 
by exposure to environmental agents such as UV and nutri- 
ent starvation or endogenously as by-products of metabolic 
pathways in different cellular compartments (e.g., respira- 
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tion in mitochondria and /^-oxidation of fatty acids in per- 
oxisomes) (Apel and Hirt 2004). To cope with increased 
ROS levels, cells evolved various means to detoxify ROS 
including antioxidant enzymes (e.g., catalases, peroxidases, 
and superoxide dismutase), and glutathione and thioredox- 
in systems (Apel and Hirt 2004; Gessler et al. 2007). Of 
these, two major antioxidant enzymes, superoxide dismu- 
tase (SOD) and catalase, play important roles in the ROS- 
scavenging mechanism. SOD catalyzes dismutation of 
superoxide to H2O2 and catalase detoxifies H2O2 to H2O 
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and O2. As part of the defense mechanism of cells, ROS was 
proposed to act as a second messenger to activate antioxi- 
dant gene transcription in response to oxidative stress and 
to regulate development and differentiation in molds, bud- 
ding, and fission yeasts, as well as plants (Toone and Jones 
1998; Apel and Hirt 2004; Miskei et al. 2009). 

In budding and fission yeasts, histidine kinases in 
two-component signaling systems function as sensors of 
oxidative stress and activate mitogen-activated protein 
kinase (MAPK) -signaling cascades, resulting in induction 
of defense gene expression (Bahn et al. 2007). In the fis- 
sion yeast Schizosaccharoniyces pombe, the basic leucine 
zipper (bZIP) transcription factor Atfl, a member of the 
ATF/cAMP response element binding (CREB) protein 
DNA-binding family, is activated by the stress-activated 
MAPK Styl/Spcl through a MAPK cascade. Atfl then pos- 
itively regulates entry of the cell into stationary growth 
phase (Moye-Rowley 2002). Expression of the cytosolic 
catalase T gene, cttl, is strongly dependent on Atfl and 
another bZIP-type transcription factor, Papl, in S. pombe 
in response to oxidative stress (Nguyen et al. 2000; Moye- 
Rowley 2002). In the budding yeast Saccharotnyces 
cerevisiae, Skol (a homolog of Atfl in S. pombe) binds to 
a cAMP-response element consensus sequence (5'-TG/ 
TACGTC/AA-3') in promoters of stress-response genes 
(Proft et al. 2005). 

Filamentous fungi use similar stress-response signaling 
pathways to those in yeast (Kawasaki et al. 2002; Furuka- 
wa et al. 2005; Hagiwara et al. 2007; Lara-Rojas et al. 
2011). Homologs of almost all components of the stress- 
activated protein kinase (SAPK) -signaling cascade in yeast 
were also found in Aspergillus nidulans (Miskei et al. 
2009). For example, AtfA, a homolog of Atfl in S. pombe, 
plays an important role in tolerance of conidia to oxida- 
tive stress in A. nidulans whUe Nap A, a homolog of Papl 
in S. pombe, is responsible for hyphal tolerance to oxida- 
tive stress (Asano et al. 2007; Hagiwara et al. 2008). In 
Aspergillus oryzae, atfA and atfB were identified and atfA 
was expressed throughout growth while atfB was 
expressed late in the growth phase and mediated oxidative 
stress tolerance in conidia generated on a solid growth 
medium (Sakamoto et al. 2008). The authors showed that 
AtfB positively regulates target genes including conidia- 
specific catA, and these genes carry CRE motifs 
(5'-TGACGTCA-3') in their promoters. Specific AtfB 
binding to these promoters was not demonstrated. 

Aspergillus parasiticus and Aspergillus flavus are two 
major producers of aflatoxin which is closely associated 
with human liver cancer. Others reported a correlation 
between ROS formation, aflatoxin production, and anti- 
oxidant enzyme activation (Fanelli et al. 1984; Narasaiah 
et al. 2006; Reverberi et al. 2006). Aflatoxin production 
in A. parasiticus was increased by oxidative stress (Jayash- 



ree and Subramanyam 2000), and it was reduced by 
antioxidants including eugenol, by stimulation of antioxi- 
dant enzyme activities such as catalase and SOD, or by 
ROS scavengers including ethylene (Jayashree and Subr- 
amanyam 1999; Reverberi et al. 2005; Huang et al. 2009). 
Surprisingly, a toxigenic strain of A. parasiticus showed 
increased activities of antioxidant enzymes such as SOD 
compared with the nontoxigenic strain and disruption of 
an SOD gene in A. flavus resulted in decreased aflatoxin 
production (Narasaiah et al. 2006; He et al. 2007). In 
A. parasiticus, disruption of ApyapA, an ortholog of yapl 
in S. cerevisiae, resulted in precocious and increased ROS 
formation and precocious aflatoxin production (Reverberi 
et al. 2007, 2008). 

Taken together, the evidence strongly suggests that afla- 
toxin biosynthesis is part of the fungal response to oxida- 
tive stress. The mechanisms that activate transcription of 
genes involved in aflatoxin biosynthesis (for clarity, we 
will designate these as aflatoxin biosynthetic genes or afla- 
toxin genes) in response to oxidative stress are not fully 
understood. However, our recent studies shed light on 
the regulatory mechanisms that integrate secondary 
metabolism and cellular response to oxidative stress (Roze 
et al. 2004, 2011). In these previous studies, we demon- 
strated that AtfB binds to promoters of seven aflatoxin 
genes which carry CRE motifs using chromatin immuno- 
precipitation (ChIP) in A. parasiticus grown under 
aflatoxin-inducing conditions (Roze et al. 2011). Electro- 
phoretic mobility shift assays (EMSA) also showed that 
AtfB binds to CREl (5'-TGACATAA-3') and AP-1 
(5'-TGAGTAC-3') binding sites in the promoter of nor-1, 
an early aflatoxin gene. 

In this study, we investigated the mechanism of 
transcriptional co-regulation of antioxidant genes 
(mycelia-specific catl and mitochondria-specific Mn sod) 
and secondary metabolism genes (the aflatoxin biosyn- 
thetic genes fas-1, ver-1, and omtA representing early, 
middle, and late aflatoxin biosynthetic genes, respectively) 
as a mean of cellular response to oxidative stress in 
A. parasiticus. We demonstrate using EMSA that the 
stress-related transcription factor AtfB binds to the pro- 
moters of early, middle, and late aflatoxin biosynthetic 
genes as well as stress-response genes that carry at least 
one CRE motif under aflatoxin-inducing conditions, but 
AtfB does not bind to the vbs promoter, which lacks a 
CRE site under the same growth condition. Based on 
nucleotide sequence and EMSA analyses, we also propose 
that three other transcription factors, MsnA, SrrA, and 
AP-1, participate in transcriptional co-regulation of 
stress-response and aflatoxin biosynthetic genes. Further- 
more, we found a conserved motif 5'-T/GNT/CA- 
AGCCNNG/AA/GC/ANT/C-3' in the promoters of fas-1, 
ver-1, nor-1, mycelial catl, and Mn sod, which extends the 
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conserved 8-base motif 5'-AGCCG/CTG/CA/G-3' identi- 
fied in the promoter of the nor-1 in our previous study; 
the core consensus sequence 5'-AAGCC-3' in the con- 
served motif is a sequence-specific-binding motif for the 
transcription factor SrrA. We also found a motif 
(5'-CCCCT-3') for a stress-response element (STRE) in 
the promoters of the aflatoxin biosynthetic genes, fas-1 
and ver-1, as well as antioxidant genes, mycelial catl and 
Mn sod. Competition EMSA analysis suggests that MsnA, 
which binds to STREs, is one player in transcriptional 
co-regulation. These data contribute to our understanding 
of the molecular mechanisms of transcriptional activation 
of aflatoxin biosynthetic genes and stress-response genes 
in response to oxidative stress in A. parasiticus. 

Experimental Procedures 

Strains, growth media, and growth 
conditions 

Aspergillus parasiticus SU-1 (ATCC 56775) was used as a 
wild-type strain for aflatoxin production in this study. YES 
liquid medium (2% yeast extract, 6% sucrose; pH 5.8) was 
used as an aflatoxin-inducing medium. A 100 mL of YES 
medium in a 250-mL flask with five 6-mm glass beads was 
inoculated with 1 x 10* spores mL^' and incubated at 30° 
C in the dark with shaking at 150 rpm as described previ- 
ously (Roze et al. 2011). 

Generation of anti-AP-1 antibodies 

Design of peptide antigens 

A BLAST search of the A. flavus genome database (http:// 
www.aspergillusflavus.org/genomics) using a partial 
sequence of A. parasiticus AP-1 (Reverberi et al. 2008) 
identified a single gene AFLA_1 29340 corresponding to a 
putative bZIP transcription factor, AP-1. Peptide antigens 
KQQNDKQRSQAKTDPD (amino acids 34-49; peptide 
102) and VNQKDVEDIMGRVK (amino acids 571-584; 
peptide 103) were designed using Antigen Profiler™ 
(Thermo Fischer Scientific, Open Biosystems, Huntsville, 
AL). 

Immunization of rabbits 

Synthetic peptides 102 and 103 carrying Cys on the C 
terminus were conjugated to the carrier protein keyhole 
limpet hemocyanin and injected into New Zealand White 
rabbits as described previously (Roze et al. 2011). 
Specificity and titer of antisera were determined for each 
peptide using ELISA. 
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Immunoaffinity antibody purification 

Polyclonal antibodies were purified from crude antiserum 
using an affinity column containing synthetic peptide 
coupled to Sepharose beads through an N-terminal 
cystein residue. Antibodies bound to the column were 
eluted into neutralizing 1 mol/L Tris buffer using a step- 
wise pH gradient. Purified antibodies were stored in anti- 
microbial saline (pH 8.2) at -H4°C or -80°C. 

Western blot analysis 

Aspergillus parasiticus SU-1 was grown in YES medium 
for designated periods of time and mycelia were harvested 
by filtration through Miracloth. Mycelia were frozen in 
liquid nitrogen. Frozen mycelia were ground in liquid 
nitrogen with a mortar and pestle, and resuspended at a 
1:1 (w/v) ratio in TSA buffer (0.01 mol/L Tris, 0.15 mol/ 
L NaCl, 0.05% NaNj, pH 8.0) containing one tablet of 
Complete mini protease inhibitor mixture (Roche Diag- 
nostics, Indianapolis, IN), 50 /(L of proteinase inhibitor 
mixture (Sigma-Aldrich, St. Louis, MO), and 0.125 
mmol/L phenylmethylsulfonyl fluoride per 10 mL. Total 
protein extracts (90 /ig per lane) or enriched nuclear pro- 
tein extracts (200 /.ig per lane) were separated by electro- 
phoresis on 12% SDS-polyacrylamide gels, transferred 
onto polyvinylidene difluoride (PVDF) membranes (Perk- 
inElmer Life Sciences, Waltham, MA), and probed with 
anti-AP-1 antibodies (102 or 103; 5 jig mL^') as a pri- 
mary antibody. The membranes were incubated with goat 
anti-rabbit secondary antibody conjugated to the fluores- 
cent tag IRDye 680 (Li-Cor Biosciences, Lincoln, NE). 
Protein bands were detected using an Odyssey infrared 
imaging system (Li-Cor Biosciences) at 795 nm. 

Preparation of cell extracts enriched in 
nuclear proteins for electrophoretic 
mobility shift assays 

Preparation of cell extracts enriched for nuclear proteins 
from A. parasiticus SU-1 for EMSA was performed as 
described previously (Roze et al. 2004). Briefly, 100 mL 
of YES medium was inoculated with conidiospores 
(1 X lO" mL"') and incubated for 24, 48, or 60 h at 30°C 
with shaking at 150 rpm. The mycelia were harvested 
by filtration through Miracloth, ground in liquid nitrogen 
using a mortar and pestle, and resuspended in lysis buffer 
as described previously. The proteins were precipitated by 
ammonium sulfate (10% and then 70%), pelleted by cen- 
trifugation, resuspended in dialysis buffer, and dialyzed 
against the dialysis buffer. The dialyzed protein solution 
was aliquated and stored at — 80°C. 
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Electrophoretic mobility shift assays 

Electrophoretic mobility shift assays, competition EMSA, 
and shift inhibition EMSA were performed as described 
previously (Roze et al. 2004). Double-stranded DNA 
fragments derived from the promoter region of /as- J, ver- 
1, omtA, vbs, mycelial catl, and Mn sod were generated 
by PGR and gel purified using a Wizard SV gel and PGR 
cleanup system (Promega, Madison, WI). Primers used to 
generate PGR fragments are shown in Table SI. The DNA 
fragments were 5' end-labeled with [T-^^P]ATP usme 



USB OptiKinase (Affymetrix, Gleveland, OH), purified 
using a Micro Bio-Spin P-30 chromatography column 
(Bio-rad, Hercules, GA), and used as probes. Two alternative 
methods were used to generate double-stranded DNA 
competitor fragments: (1) A 51 bp NorR4 fragment was 
generated by PGR using pUG19GUS as a template and gel 
purified using a Wizard SV gel and PGR cleanup system 
(Promega) (Roze et al. 2011). Primers used to generate 
NorR4 are shown in Table SI. (2) A 28-bp Fas fragment 
including the motif 5'-AGGGGTGA-3' (Roze et al. 2011) 
derived from the fas2/fasl intergenic region was generated 



(A) 



(B) 



CRE3 CRE2 STRE2 



CREl STREl 



fas-2lfas-I intergenic region (701 bp) 



CRE3 CRE2 STRE2 



-i- 



Fas2 (298 bp) 



CREl STREl 



Fasl (399 bp) 



Fasl-b(230 bp) 

CREl STREl 



Fasl-a(189bp) 



(C) 



CRE5 CRE4CRE3 



STREl CRE2 CREl 



(E) 



CRE3 STRE2 



Fas (28 bp) 



STRE1CRE2CRE1 



\ 



1 1 mcatl 

mycelial catl promoter region (Mcatl, 577 bp) upstream ATG 
CRE3 STRE2 



Mcatl-b(339 bp) 



STRE1CRE2CRE1 



norAlver-l intergenic region (81 1 bp) 
CRE5 CRE4CRE3 



Verl-b(311 bp) 



STREl CRE2 CREl 



(D) 



Verl-a (499 bp) 
STREl CRE2 

Verl (55 bp) 
STREl 
-I- 
Verlm (55 bp) 



CRE9CRE8 CRE7CRE5 STRE 1CRE5 CRE4CRE3CRE2CRE 1 
Ill 1 1 I I I I 



omtAlordA intergenic region (1320 bp) 
CRE9CRE8CRE7 



OintA (596 bp) CRE5 STRE1CRE5 CRE4CRE3CRE2CRE1 



OrdA (692 bp) 



Mcatl-a(259bp) 



(F) 



CRE3 STRE2 CRE2CRE1 STREl 



I \ 



Mn sod promoter region (Msoti, 455 bp) upstream ATG 



Figure 1. Promoter regions of aflatoxin biosynthetic and antioxidant genes used in electrophoretic mobility shift assays. (A) Schematic 
representation of the relevant genes in the aflatoxin gene cluster. (B) Schematic of fas-2/fas-1 intergenic region. The fas-2/fas-1 intergenic region 
was divided into two fragments designated Fas2 and Fasl . The Fasl fragment was further subdivided into two smaller fragments designated 
Fasl-b and Fasl-a. Five putative c/s-acting sites are shown including CRE and stress-response element (STRE) sites. The position of a 28-bp Fas 
competitor is shown. (C) Schematic of norA/ver-1 intergenic region. The norA/ver-1 intergenic region was divided into two fragments designated 
Verl-b and Verl-a. Six putative c/s-acting sites are shown including CRE and STRE sites. The position of 5B-bp Verl and Verl m competitors is 
shown. (D) Schematic of omtAlordA intergenic region. The omtAJordA intergenic region was divided into two fragments designated OmtA and 
OrdA. Ten putative c/s-acting sites are shown including CRE and STRE sites. (E) Schematic of mycelial caf7 promoter region designated Mcatl. 
The mycelial cat; promoter region was divided into two fragments designated Mcatl -b and Mcatl -a. Five putative c/s-acting sites are shown 
including CRE and STRE sites. (F) Schematic of Mn sod promoter region designated Msod. Five putative c/s-acting sites are shown including CRE 
and STRE sites. 
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by annealing of two complementary single-stranded syn- 
thetic oligonucleotides (Fig. IB). The sequences used to 
generate the Fas fragment are shown in Table SI. This same 
strategy was used to generate a 55-bp Verl fragment, which 
carried two motifs (STREl: CCCCT and CRE2: TGACC- 
CAG), and a 55-bp Verl m fragment including the STRE 
motif only (Verl with CRE2 changed from TGAC to 
CTTC) (Fig. IC) [Correction added on 04 February 2013, 
after first online publication: "CREl" has been changed to 
"CRE2".]. The sequences used to generate the Verl and 
Verl m fragments are shown in Table SI. 

Electrophoresis on 5% nondenaturing polyacrylamide 
(80:1 acrylamide/bisacrylamide) gels was used to separate 
DNA-protein complexes. Twenty femtomole of each 
'^P-labeled probe was incubated for 15 min at 30°C with 
2 iig of poly(dl-dC), 7.5 fig of bovine serum albumin, 
and competitor (if desired) with 5 /ig of nuclear protein 
extract (added last) in a DNA-binding buffer (15% glyc- 
erol, 15 mmol/L HEPES [pH 7.9], 100 mmol/L KCl, 
1 mmol/L EDTA, 2 mmol/L DTT) in 25 fiL of final bind- 
ing reaction. After separation, the gels were dried and 
exposed to x-ray film. 

Polyclonal anti-AtfB antibodies (YSR) were used to 
block formation of DNA-protein complexes (shift inhibi- 
tion EMSA) as described previously (Roze et al. 2011). 
Five micrograms of nuclear protein extracts were incu- 
bated in a DNA-binding buffer for 15 min at 30°C with 
5 of anti-AtfB antibody (5.5 ^ig /^(L^') or preimmune 
serum obtained from the same rabbit as the antibody was 
generated. Then, a ^^P-labeled probe was added and the 
mixture was incubated for an additional 15 min at 30°C. 
Finally, DNA-protein complexes were resolved by electro- 
phoresis. 

In some reactions, rabbit polyclonal antibodies specific 
for AP-1 were used to block formation of DNA-protein 
complexes (shift inhibition EMSA) as described above. 

DNA sequence analysis of promoter regions 
in A. parasiticus 

DNA sequence analysis of PGR products carrying pro- 
moter regions of fas-1, ver-1, otntA, vbs, mycelial catl, 
and Mn sod was conducted at the Research Technology 
Support Facility (RTSF) at Michigan State University. 
Comparison of promoter sequences in A. parasiticus and 
A. flavus was performed using EMBOSS pairwise 
sequence alignment. 

Total RNA isolation and transcript level 
analysis by quantitative real-time PCR 

Transcript levels in fungal cells were analyzed using real- 
time PCR as described previously (Roze et al. 2011). 



Aspergillus parasiticus SU-1 was cultured in YES liquid 
medium, and the mycelia were harvested and frozen in 
liquid N2 at 24, 48, or 60 h. Total RNA was isolated 
from two biological replicates by the TRIzol method 
(TRIzol Reagent; Invitrogen, Carlsbad, CA). Total RNA 
was treated with RNase-free DNase (Qiagen, Valencia, 
CA) and RNA quality was examined by an Agilent 2100 
Bioanalyzer (Agilent Technologies, Santa Clara, CA). For 
reverse transcription, 1 /ig of total RNA was treated with 
gDNA WipeOut, and cDNA was prepared with the 
QuantiTect reverse transcription kit (Qiagen). Primers 
were designed using Primer Express software 3.0 
(Applied Biosystems, Carlsbad, CA) and PrimerQuest^'^ 
(Integrated DNA Technologies, CoralviUe, lA). Primers 
used for quantitative real-time PCR (qRT-PGR) analysis 
are listed in Table S2. For qRT-PGR analysis, 25 ng of 
reverse-transcribed RNA was amplified using an ABI 
PRISM 7900HT sequence detection system with Power 
SYBR Green PGR Mastermbc (Applied Biosystems). 
Reactions were carried out using an initial incubation of 
2 min at 50°G and 10 min at 95°C. This was followed 
by 40 cycles consisting of 15 sec at 95°C (denaturation) 
and 1 min at 60°C (annealing and extension). Resulting 
data were analyzed by the relative quantification method 
using SDS 2.1 software (Applied Biosystems). Relative 
gene expression levels were determined by the standard 
curve method (Applied Biosystems). The standard curve 
was derived by plotting cycles to threshold values versus 
the logarithm of known concentration of reverse- 
transcribed RNA between the ranges of 0.001 and 100 ng 
per reaction. For each gene, the relative level of mRNA 
at each time point is represented by the fold change of 
mRNA levels of a target gene divided by mRNA levels 
of /J-tubulin at the same time point. Statistical analysis 
for all qRT-PGR was performed using SigmaPlot® soft- 
ware (v 11.0; Systat Software Inc., San Jose, CA). qRT- 
PGR data are presented as mean ± SE, n = 4, and were 
analyzed by one-way repeated measure analysis of 
variance. 

Results 

Preliminary EMSA studies detected DNA-protein com- 
plex formation between EMSA protein extracts from 
A. parasiticus grown in YES aflatoxin-inducing medium 
for 48 h and DNA fragments derived from the intergen- 
ic or promoter regions of the divergently transcribed 
fas-2/fas-l genes (the entire intergenic region, 701 bp), 
the norAlver-1 genes (the entire intergenic region, 
811 bp), the divergently transcribed omtA/ordA genes 
(the entire intergenic region, 1320 bp), the mycelial catl 
gene (577 bp upstream from ATG), and the Mn sod 
gene (455 bp upstream from ATG) (Fig. 1). In order to 
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further characterize these DNA-protein interactions, 
these large DNA fragments were subdivided into smaller 
pieces. The fas-2/fas-l intergenic region was divided into 
Fas2 (298 bp) and Fasl (399 bp). Fasl was further split 
into Fasl-b (230 bp) and Fas I -a (189 bp). The norAI 
ver-1 intergenic region was divided into Verl-b (311 
bp) and Verl-a (499 bp). The omtAlordA intergenic 
region was divided into OmtA (596 bp) and OrdA 
(692 bp). The Mcatl fragment (577 bp upstream from 
ATG) was divided into Mcatl-b (339 bp) and Mcatl-a 
(259 bp) (Fig. 1). Fragments Fasl-a, Verl-b, OmtA, 
Mcatl-a, and Msod demonstrated prominent DNA- 
protein complexes in EMSA, and therefore, we contin- 
ued to employ those fragments in EMSA experiments 
described below. 

Time-course of DNA-protein complex 
formation in aflatoxin biosynthetic and 
stress-response gene promoters 

To correlate transcription factor-binding activity with 
activation of gene expression, EMSA was conducted on 
Fasl-a, Verl-b, OmtA, Mcatl-a, and Msod promoter 
fragments using nuclear protein extracts from A. parasiti- 
cus grown in YES aflatoxin-inducing medium for 24 h 
(prior to aflatoxin gene expression), 48 h (during maxi- 
mal rates of aflatoxin gene expression), or 60 h (after 
maximal rates of aflatoxin gene expression) (Trail et al. 
1995). DNA-protein complexes were detected with Fasl- 
a, Verl-b, and OmtA promoter fragments using the 24-h 



(A) Fasl-a Verl-b OmtA (B) Mcatl-a Msod 




SU-1 (5(ig) -+ + + -+ + + -+ + + SU-1 (5 ^g) - + + + - + + + 
Time(h) 244860 244860 244860 Time (h) 244860 244860 

Figure 2. Time-course of DNA-protein complex formation as 
detected by electropiioretic mobility shift assays using aflatoxin (fas- 1 
and ver-1) and antioxidant gene (mycelial cat? and Mn socO 
promoters. Aspergillus parasiticus SU-1 was grown for 24, 48, or 
60 h at 30°C in YES medium. Cell extracts enriched in nuclear 
proteins were prepared as described in Experimental Procedures. Five 
micrograms of enriched nuclear protein extracts was added to a 
P^^-labeled promoter probes for each aflatoxin or antioxidant gene. 
(A) Fasl-a, Verl-b, and OmtA probes. (B) Mcatl-a and Msod probes. 



nuclear protein extract; similar sized but more abundant 
DNA-protein complexes and new complexes were 
observed with the same promoter fragments using 
48- and 60-h nuclear protein extracts (Fig. 2A). These 
data are consistent with strong activation of the promot- 
ers of the aflatoxin biosynthetic genes at 48 and 60 h 
observed under aflatoxin-inducing conditions (Roze et al. 
2011). The promoter of the stress-response gene, Mn sod, 
showed a similar pattern of DNA-protein complex for- 
mation as the aflatoxin biosynthetic genes (Fig. 2B). In 
contrast, for the promoter of the stress-response gene 
Mcflfi (Mcatl-a), the most abundant DNA-protein com- 
plexes were observed using the 24-h nuclear protein 
extract while DNA-protein complexes were observed at 
significantly reduced levels and some complexes 
disappeared using 48- or 60-h nuclear protein extracts 
(Fig. 2B). These data indicate that Mn sod and the afla- 
toxin genes /fls-i, ver-1, and omtA are co-regulated while 
Mcflfi is differentially expressed early in growth. 

Identification of c/s-acting-binding sites in 
the early and middle aflatoxin biosynthetic 
and stress-response gene promoters 

To locate specific cis-acting sites in promoters of fas-1, ver- 
1, mycelial catl, and Mn sod, competition EMSA was per- 
formed. We previously demonstrated that a 51 -bp NorR4 
fragment derived from the nor-1 promoter, contained CRE 
and AP-1 sites (5'-TGACATAA-3' for CRE; and 5'-TGAG- 
TAC-3' for AP-1) and it competed efficiently for complex 
formation with a 170-bp NorR promoter fragment using a 
48-h nuclear protein extract (Roze et al. 2011). We 
hypothesized that the promoters of other aflatoxin biosyn- 
thetic genes including fas-1 and ver-1 and the stress- 
response genes mycelial catl and Mn sod are activated in 
response to oxidative stress in a similar way by binding 
identical transcription factors to analogous cfs-acting ele- 
ments in the promoters. Thus, first we used the 51 -bp 
NorR4 fragment as a competitor in EMSA. The 51 -bp 
NorR4 fragment containing both sites competed efficiently 
for DNA-protein complexes formed in Fasl-a, Verl-b, 
Mcatl-a, and Msod promoter fragments using a 48-h 
nuclear protein extract (Fig. 3A and B). These data suggest 
that promoters of early (fas-1) and middle (ver-1) aflatoxin 
genes as well as stress-response genes (mycelial catl, and 
Mn sod) carry similar CRE sites, AP-1 sites, or both sites as 
observed in the nor-1 promoter. 

We previously identified a conserved 8 -base motif 
(5'-AGCCG/CTG/CA/G-3') within promoters of eight 
aflatoxin genes by sequence-based MEME (a computer 
software package for sequence-based motif analysis of 
nucleotide sequences) motif analysis and showed AtfB 
binding to five of these eight promoters including fas-1 
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(A) Fasl-a Verl-b (B) Mcatl-a Msod (C) Fasl-b Fasl-a Verl-b 




su-1 (5 ^g) - + + + -+ + + su-i (5 ng) - + + + - + + + su-1 (5^g) - + + + _ + + + _ + + + 

Time(h) 48 48 Time (h) 48 48 Time (h) 48 48 48 



Figure 3. Competition electrophoretic mobility shift assays. Aspergillus parasiticus SU-1 was grown for 48 h at 30°C in YES medium. (A) Fasl-a 
and Verl-b probes. (B) Mcatl-a and Msod probes. The Bl-bp nonlabeled NorR4 fragment (contains CRE and AP-1 sites; 50- and 250-fold molar 
excess) was added to compete for labeled Fasl-a, Verl-b, Mcatl-a, or Msod probes. (C) Fasl-b, Fasl-a, and Verl-b probes. The 28-bp 
nonlabeled Fas fragment (contains 8-base motif, 5'-AGCCG/CTG/CA/G-3'; 50- and 250-fold molar excess) was added to compete for labeled 
Fasl-b, Fasl-a, or Verl-b probes. 



and nor-1 using ChIP analyses; the data suggested that 
the conserved 8-base motif may contribute to AtfB bind- 
ing (Roze et al. 2011). Therefore, we tested whether the 
conserved 8-base motif actually binds AtfB or another 
transcription factor within the fas-1 (which contains the 
motif) and ver-1 (which does not contain the motif) 
promoters using competition EMSA. A 28-bp Fas frag- 
ment containing the 8-base motif, which was derived 
from the Fasl-b promoter fragment (Fig. IB), competed 
weakly for DNA-protein complex formation in Fasl-b, 
Fasl-a, and Verl-b promoter fragments using 48 -h 
nuclear protein extract (Fig. 3C). The data suggest that 
the ver-1 promoter carries a motif similar to the con- 
served 8-base motif in the fas-1 promoter and that AtfB 
and/or another transcription factor binds specifically to 
this motif in the promoters of fas-1 and ver-1. 

AtfB binds in aflatoxin biosynthetic and 
stress-response gene promoters that 
carry a CRE motif 

Previously, we demonstrated that AtfB binds to the nor-1 
promoter using EMSA protein extracts from A. parasiticus 
grown for 48 h in YES aflatoxin-inducing medium and 
antibodies that specifically recognize AtfB (Roze et al. 
2011). In this study, first, we determined whether AtfB 
binds early, middle, and late aflatoxin gene promoters 
that carry CRE sites (see Table 1) using nuclear protein 
extracts prepared from A. parasiticus grown in YES med- 
ium for 24, 48, or 60 h and AtfB-specific polyclonal anti- 
bodies to conduct shift inhibition EMSA. 



Table 1. CRE sites in aflatoxin and antioxidant gene promoters used 
in this study. 



Aflatoxin and 


Total number of all 




antioxidant gene 


CRE sites 


All CRE sites containing 


promoters 


containing TGAC 


TGAC^ 


fas-21 fas-1 


3 


TGACGGAC, TGACCGAT, 


(701 bp) 




TGACAAGC 


norAlver- 1 


5 


TGACGAGA, TGACTGAG, 


(811 bp) 




TGACCAAG, 






TGACCCAG, TGACGGGA 


omtA/ordA 


9 


TGACTTCG, TGACGCAA, 


(1320 bp) 




TGACTCGT, TGACCATA, 






TGACGCCA, 






TGACGAAC, 






TGACCCAC, TGACGCCC, 






TGACGATC 


ordA/vbs 


None 


None 


(511 bp) 






Mycelial cat? 


3 


TGACCTAC, TGACACTT, 


577 bp 




TGACAGTT 


upstream ATG 






Mn sod 455 bp 


3 


TGACTGAC, TGACTATG, 


upstream ATG 




TGACCAGA 



'The nucleotides showing the same sequence with the consensus CRE 
motif (TGACGTCA) were indicated in bold. 



DNA-protein complexes were observed with Fasl-a, 
Verl-b, and OmtA promoter fragments using 24-, 48-, or 
60-h nuclear protein extracts and complex formation was 
inhibited by pretreatment of the nuclear protein extracts 
with anti-AtfB antibodies before addition of the promoter 
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SU-l (5 fig) - + + + - + + + - + + + SU-1 (5|ig) - + + + - + + + - + ++ Sy-I (5 fig) - + + + 

Time(h) 24 48 60 Time (h) 24 48 60 Time (h) 48 




Figure 4. Electrophoretic mobility sinift assays of AtfB binding in the fas-1, ver-1, omtA, vbs, mycelial cat?, and Mn sod promoters. Aspergillus 
parasiticus SU-1 was grown for 24, 48, or 60 h at 30°C in YES medium. Enriched nuclear protein extracts were prepared as described in 
Experimental Procedures. Five micrograms of enriched nuclear protein extracts was added to a P^^-labeled promoter probe for each aflatoxin or 
antioxidant gene. Anti-AtfB antibodies (VSR) or preimmune serum was added to determine whether these could block protein/DNA interaction 
(shift inhibition). (A) Fasl-a probe. Nuclear protein extracts were used from 24-, 48-, or 60-h culture. (B) Verl-b probe. Nuclear protein extracts 
were used from 24-, 48-, or 60-h culture. (C) OmtA probe. Nuclear protein extracts were used from 48-h culture. (D) Vbs probe. Nuclear protein 
extracts were used from 48-h culture. (E) Mcatl-a probe. Nuclear protein extracts were used from 24-, 48-, or 60-h culture. (F) Msod probe. 
Nuclear protein extracts were used from 24-, 48-, or 60-h culture. 



fragment but not by preimmune serum obtained from the 
same rabbit that the antibody was generated (Fig. 4A-D). 
These data suggest that complex formation is due to spe- 
cific AtfB binding to the promoter fragments at all time 
points analyzed. The Verl-b probe produced two DNA- 
protein complexes (upper and lower) as shown in 
Figures 2A, 3C, and 4B. Only the upper complex contains 
AtfB, because the complex was competed by anti-AtfB 
(Fig. 4B); the nature of the protein(s) in the lower band 
is unknown. Similarly, the Fasl-b promoter fragment 
demonstrates two complexes as shown in Figure 3C. Only 
the upper complex with higher intensity competed by 
Fas-unlabeled probe. These findings are also consistent 
with previous results obtained with anti-AtfB antibodies 
in Chip analysis that showed enrichment of AtfB binding 
to each promoter carrying CRE sites under aflatoxin- 
inducing conditions (Roze et al. 2011). 



As a control, we used a vbs promoter fragment (Vbs: 
508 bp), which lacks CRE sites (Table 1) (Roze et al. 
2004), for EMSA. Vbs-protein complex formation was 
not observed using a 48-h nuclear protein extract 
(Fig. 4D). These data confirm previous results obtained 
using Chip analyses, in which the vhs promoter did not 
bind AtfB at any time point under aflatoxin-inducing 
conditions (Roze et al. 2011). 

When taken together, these results indicate that AtfB 
binds at early (fas-1), middle (ver-1), and late (omtA) 
aflatoxin biosynthetic gene promoters that carry CRE 
sites. 

Next, to determine whether AtfB also binds at stress- 
response gene promoters that carry CRE sites (see 
Table 1), shift inhibition EMSA was performed on the 
mycelial catl and Mn sod promoters. A Mcatl-a-protein 
complex was observed in EMSA using 24-, 48-, or 60-h 
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nuclear protein extracts (Fig. 4E) and complex formation 
was inhibited by pretreatment of the nuclear protein 
extracts with anti-AtfB antibodies similar to results 
observed for the aflatoxin biosynthetic genes (Fig. 4E). 
For a second stress-response gene, an Mn sod promoter 
(Msod, 455 bp)-protein complex was detected by EMSA 
using a 24-, 48-, or 60-h nuclear protein extracts (Fig. 4F) 
and complex formation was inhibited by pretreatment of 
the nuclear protein extracts with anti-AtfB antibody using 
48- or 60-h nuclear protein extracts. However, antibody 
pretreatment did not prevent complex formation when 
the 24-h protein extract was used (Fig. 4F). These results 
imply that the DNA-protein complex formation on the 
Mn sod promoter depends on AtfB at 48 and 60 h, but 
not at 24 h. 

AP-1 -specific polyclonal antibodies do not 
block DNA-protein complex formation 

As previous data suggested that an AP-1 site (5'-TGAG- 
TAC-3') played a role in activation of the nor-1 promoter, 
we tested whether DNA-protein complex formation in 
other aflatoxin genes and stress-response genes is inhib- 
ited by pretreatment of the nuclear protein extract with 
AP-1 -specific antibodies (shift inhibition EMSA). To 
demonstrate antibody specificity, we conducted Western 
blot analysis of the A. parasiticus nuclear protein extract. 
AP-1 antibodies detected a 64-kDa protein in the A. para- 
siticus nuclear protein extract (Fig. SI), and this mass 
agreed with the predicted molecular mass of AP-1 
(63.8 kDa) in A. flavus (GenBank accession no: 
EED48711.1). AP-1 antibody pretreatment of the nuclear 



protein extracts did not affect DNA-protein complex for- 
mation within the Fasl-a, Ver-lb, Mcatl-a, or Msod pro- 
moter fragments (Fig. S2) suggesting that AP-1 does not 
bind directly to promoters of these target genes. 

Aflatoxin biosynthetic and stress-response 
gene promoters contain a conserved 
AAGCC motif 

Because shift inhibition EMSA analysis demonstrated that 
AtfB binds specifically to Fasl, Verl, Mcatl, and Msod 
promoter fragments and that the NorR4 fragment com- 
peted efficiently for DNA-protein complexes formed in 
the promoter fragments, we applied sequence-based 
MEME motif analysis to search for motifs that were con- 
served in these promoters and in the NorR4 competitor 
fragment. A 5-base conserved motif 5'-AAGCC-3' was 
found in all five promoter regions (Fig. 5A and B), and 
this motif is a part of the 8-base motif 5'-AGCCG/CTG/ 
CA/G-3' identified in the nor-1 promoter in previous 
work (Roze et al. 2011). Sequence-based MEME analysis 
of these same five promoters extended the 8-base motif 
to a 15-base conserved motif 5'-T/GNT/CAAGCCNNG/ 
AA/GC/ANT/C-3' (Fig. 5C and D). Interestingly, the con- 
served core consensus sequence 5'-AAGCC-3' has been 
shown previously to be recognized by the transcription 
factor SrrA, an A. parasiticus ortholog of S. cerevisiae 
Skn7. Skn7 recruits Yapl to the promoters of yeast target 
genes and activates cooperatively their expression in 
response to oxidative stress (He and Fassler 2005; Mul- 
ford and Fassler 2011); transcriptional activation occurs 
through sequence-specific interaction of Skn7 with one of 



(A) 

Name Strand Start p-value Sites 

Msod + 293 2.52e-05 GCATCCATTC CCAGGCTT ATATTCTCCC 

Mcatl-a - 69 3.85e-05 AAGAGAAGTT CCTGGCTT GCAGGTACAC 

NorR4 - 26 6.52e-05 TTCATTATGT CACGGCTT AGAAAAGTAC 

Fasl-a + 23 6.52e-05 AGGAAGACGC CTCGGCTT GTCAAGAATG 

Verl-b - 187 l.24e-04 GCAAAGGTTC TCAGGCTT ATAATATGTT 



(C) 



Name 


Strand 


Start 


p-value 


Sites 


Msod 


+ 


3 


2.10e-09 


GAA TATAAGCCTGGGAAT 


NorR4 


+ 


3 


4.98e-07 


TTT TCTAA CC ToACAT 


Verl-a 


+ 


3 


4.98e-07 


TAT TATAA CCT -A-AAC 


Mcatl-a 


+ 


3 


8.32e-07 


ACC TGCAAGCCAGGAACT 


Fasl-a 




0 


1.18e-06 


CTT GACAAGCCGAGGCGT 




Figure 5. Sequence-based motif analysis of fas-1, ver-1, mycelial cati, and Mn sod promoters using MEfVlE. Fasl-a, Verl-b, Mcatl-a, and Msod 
promoters fragment and the NorR4 competitor used for competition electrophoretic mobility shift assays were analyzed for conserved motifs 
using sequence-based MEME motif analysis. A conserved motif (B'-T/GNT/CAAGCCNNG/AA/GC/ANT/C-3') was found in the promoter regions, 
where the core consensus sequence was B'-AAGCC-3'. (A) Sites carrying the core consensus sequence in five promoter fragments. (B) Sequence 
logo of the core consensus sequence. (C) Sites carrying the conserved motif in five promoter fragments. (D) Sequence logo of the consen/ed 
motif. 
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Table 2. AP-1 sites in afiatoxin and antioxidant gene promoters used 
in tliis study. 



Aflatoxin and 


Total number of all AP-1 




antioxidant gene 


sites witlnin 60 bp near 


All AP-1 sites within 


promoters 


SrrA site 


60 bp near SrrA site' 


Fasl-a (189 bp) 


2 


TGAGTGA, 






TGAGGCC 


\/pr1-h ^^1 1 hn^ 


4 


TTAGTCC 






TGAGCCT, 






TGAGAAC, 






TTACGTC 


Mcat1-a (259 bp) 


1 


TGAGTTC 


Msod (455 bp) 


4 


TGAGGAC, 






TGAGGGA, 






TTACAAC, 






TGAGGGA 



'The nucleotides showing the same sequence with the consensus 
AP-1 motif (TT/GAC/GT/AA/CA) were indicated in bold. 

four ris-regulatory elements within promoters (5'- 
GGCNGGC-3', 5'-GGCNNGGC-3', 5'-GGCNAGA-3', or 
5'-GGCNNAGA-3'). We observed that the AP-1 and SrrA 
motifs are located in close proximity to each other, lying 
within 60 bp in the promoters of aflatoxin {fas-1 and 
ver-1) and stress-response genes (mycelial catl and Mn 
sod) (Table 2). Our findings are consistent with the 
observation that Yapl and Skn7 recognition sites in S. ce- 
revisiae are located within 20-50 bp of each other (He 
and Fassler 2005). 



(A) 



(B) 



Competitor Competitor ^ofo-^^^^c^^^lsgc^V^ 



DNA-protein 
complexes 




DNA-protein 
complexes 




su-1 (5 ^g) - -n- + - -n- + su-i (5 ^g) - + ++-+ + + 

Time(h) 24 60 Time (h) 24 60 

(C) Mcatl-a 



(D) 



Competitor ^o^o^c^^^^^o^o^A'^ Competitor ^o^o^c^^^^^^^^ 



DNA-protein 
complexes 





DNA-protein 

complexes ^ M b *^ || | 



llllllll 



SU-l (5 fig) - + + -f - + -f + SU-1 (5 ^g) - -1- + + - + -f ■ 



Time (li) 



24 60 



Time (li) 24 60 



MsnA binds promoters of aflatoxin 
biosynthetic and stress-response genes 
that carry a STRE motif 

Msn2p and Msn4p (C2H2-type zinc-finger regulators) 
bind to a conserved STRE (5'-CCCCT-3') in promoters of 
stress-response genes (e.g., cytosolic catalase gene, CTTl) 
in response to oxidative stress (H2O2) and activate their 
expression in S. cerevisiae (Marchler et al. 1993; Marti- 
nez-Pastor et al. 1996). The Msn2p ortholog, MsnA, was 
identified in A. parasiticus and A. flavus (Chang et al. 
2011). Because of the possible link between MsnA and 
regulation of stress response, we searched for a STRE 
motif in promoters of aflatoxin and stress-response genes 
in A. parasiticus; we identified a STRE motif in the pro- 
moters of fas-1, ver-1, mycelial catl, and Mn sod (Fig. 1). 
To test whether the conserved STRE motif is involved in 
DNA-protein complex formation at these promoters, we 
conducted competition EMSA. We used a 55-bp DNA 
fragment derived from the ver-1 promoter which contains 
both STRE and CRE motifs (5'-CCCCTGGGTCA-3') 
(Fig. IC) to compete for DNA-protein complex formation. 
The 55-bp Verl fragment competed for DNA-protein 



Figure 6. Competition electrophoretic mobility shift assays of fas-7, 
ver-1, mycelial caf?, and Mn sod promoters using a 55-bp Verl 
promoter fragment as a competitor. Aspergillus parasiticus SU-1 was 
grown for 24 or 60 h at 30°C in YES medium. The 55-bp nonlabeled 
Verl promoter fragment (contains both stress-response element and 
CRE sites; 50- and 250-fold molar excess) was added to compete for 
labeled Fasl-a, Verl-b, Mcatl-a, or Msod probes. (A) Fasl-a probe. 
(B) Verl-b probe. (C) Mcatl-a probe. (D) Msod probe. 

complex formation in Fasl-a, Ver-lb, and Mcatl- 
a-promoter fragments using 24- or 60-h nuclear protein 
extracts (Fig. 6A-C). However, the 55-bp Verl promoter 
fragment did not compete for DNA-protein complex for- 
mation in the Msod promoter fragment using the 24- or 
60-h nuclear protein extracts (Fig. 6D). To identify com- 
plexes dependent on STRE only, competition EMSA was 
performed using a 55-bp Verl m fragment which carries 
the STRE motif and a mutated version of the CRE motif 
(TGAC changed to CTTC) (Fig. IC). EMSA using the 
55-bp Verl m fragment as a competitor showed the same 
results as those using the 55-bp Verl fragment (Fig. S3). 
These data suggest that a DNA-protein complex binds 
the STRE motif in fas-1, ver-1, and mycelial catl promot- 
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Verl Verlm 




SU-1 (5 fig) -+ + +-+ + + 
Time (h) 24 48 60 24 48 60 



Figure 7. Time-course of DNA-protein complex formation as 
detected by electroplnoretic mobility sh\h assays using oligonucleotide 
probes carrying a stress-response element (STRE) site. Aspergillus 
parasiticus 5U-1 was grown for 24, 48, or 60 li at 30°C in YES 
medium. Enriched nuclear protein extracts were prepared as 
described in Experimental Procedures. Five micrograms of enriched 
nuclear protein extracts was added to the P^^-labeled 55-bp Verl 
(contains STRE and CRE sites) or Verl m (contains STRE site) 
oligonucleotide probes in each lane. 

ers, and this complex contains MsnA. However, MsnA is 
not part of the DNA-protein complex that forms in the 
Mn sod promoter. 

To analyze time dependence of DNA-protein complex 
formation in the 55-bp Verl and Verl m fragments, 
EMSA was conducted using 24-, 48-, and 60-h nuclear 
protein extracts. We observed two major DNA-protein 
complexes (upper and lower) and minor DNA-protein 
complexes (Fig. 7). The intensity of the two major DNA- 
protein complexes at 24 h was maximal and decreased at 
48 and 60 h (Fig. 7). The timing of appearance and 
quantity of the lower DNA-protein complex on the 
55-bp Verl promoter fragment was similar to that of the 
Verl m fragment, supporting the idea that formation of 
the lower DNA-protein complex is mediated though the 
STRE motif However, formation of the upper DNA- 
protein complex was dependent on the CRE site. 
Formation of complexes on both Verl and Verl m pro- 
moter probes was parallel to the expression pattern of 
mycelial catl described above (Fig. 2B). Interestingly, we 
observed one minor DNA-protein complex located 
between the two major DNA-protein complexes. How- 
ever, it appears the minor complex formation was not 
affected by the mutation in the Verl m fragment. 

These data suggest that MsnA is already bound to 
promoters of stress-response genes such as mycelial catl 
at 24 h in order to activate gene expression in response 
to oxidative stress (H2O2). However, although MsnA 



binds to promoters of aflatoxin genes such as fas-1 and 
ver-1 at 24 h, it is unable to activate these genes (Figs. 6, 
8). MsnA activation of aflatoxin genes could require 
activity of additional transcription factors such as AflR 
and AtfB. 

Expression of genes encoding aflatoxin 
biosynthetic and antioxidant enzymes and 
stress-related transcription factors 

We analyzed levels of transcripts of early and middle afla- 
toxin biosynthetic genes {fas-1 and ver-1), antioxidant 
genes (mycelial catl and Mn sod), and oxidative stress- 
related transcription factor genes {atfB, AP-1, and msriA) 
in A. parasiticus SU-1 grown in aflatoxin-inducing YES 
medium at 24, 48, or 60 h using real-time PGR. Tran- 
scripts oi fas-1, ver-1, and Mn sod were barely detected at 
24 h but increased significantly at 48 h and stayed 
unchanged at 60 h (Fig. 8A, B, and D). This pattern of 
expression is consistent with the pattern of expression of 
aflatoxin genes observed previously (Roze et al. 2011). In 
contrast, the antioxidant gene encoding mycelial catl and 
genes encoding three oxidative stress-related transcription 
factors {atfB, AP-1, and msriA) were all expressed at rela- 
tively high levels at 24 h and remained unchanged at 48 
and 60 h (Fig. 8C and E-G). These data support the EMSA 
data shown in Figure 2 and indicate two distinct patterns 
of gene expression. We propose that expression of the 
aflatoxin genes {fas-1 and ver-1) and Mn sod is regulated 
by a similar mechanism, whereas expression of mycelial 
catl and the oxidative stress-related transcription factor 
genes {atfB, AP-1, and msriA) is regulated by a different 
mechanism. 

Discussion 

Our work for the first time provides strong evidence to 
suggest that at least four transcription factors AtfB, MsnA, 
SrrA, and AP-1 act together to regulate genes involved in 
stress response and secondary metabolism. The target 
genes are located on different chromosomes (based on 
the A. flavus genome). The aflatoxin gene cluster that 
carries fas-1 and ver-1 is located on chromosome 3, the 
antioxidant gene, mycelial catl, is located on chromosome 
5, and Mn sod is located on chromosome 2. Nucleotide 
sequence analysis demonstrated that the promoters of 
mycelial catl and Mn sod in A. parasiticus were highly 
conserved with each other and also showed 90% identity 
to the corresponding promoter regions in A. flavus (Fig. S4). 
In contrast, the promoters of aflatoxin genes in these 
two species exhibited significantly lower identity (50%) 
(Ehrlich et al. 2005). This observation may suggest that 
the regulation of expression of antioxidant genes is more 
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Figure 8. Expression of fas-1, ver-1, mycelial cat/, Mn sod, atfB, AP-1 , and MsnA. Aspergillus parasiticus SU-1 was grown for 24, 48, or 60 h at 
30°C in YES medium. Total RNA extraction and real-time PCR analyses were performed as described in Experimental Procedures. The relative level 
of mRNA is represented as the mRNA level of the target gene divided by the mRNA level of /J-tubulin at the same time point. Bars represents 
mean ± SE (n = 4). Statistical analysis was performed by one-way repeated measure analysis of variance. Same lowercase letter indicates that 
there is no statistically significant difference between two measurements. Different lowercase letters specify a statistically significant difference 
between two measurements. P values are as follows: fas-1 and ver-1 (P = 0.002), mycelial cafi and atfB (P = 0.011), Mn sod (P < 0.001), and 
AP-1 and MsnA (P = 0.015). Two independent biological replicates were performed showing the same trend. Two duplicates samples were 
analyzed for each biological replicate. 
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Table 3. Comparison of transcription factors wfiicfi bind to CRE and AP-1 sites in Aspergillus oryzae, Aspergillus nidulans, Saccharomyces cerevisi- 
ae, Schizosaccharomyces pombe, and mammalian cells 



Transcription factors 


A. oryzae 


A. nidulans 


5. cerevisiae 


S. pombe 


Mammalian cells 


Binding to CRE sites 


AlfA 


AtfA 


Skol(Acrl) 


Atfl 


ATF-2 


Binding to CRE sites 


AtfB 


NR^ 


NR 


NR 


ATF-1 


Binding to AP-1 sites 


NR 


NapA 


Yapl 


Papl 


c-Jun 



'NR represents that it is not reported in literature. 



highly conserved evolutionarily than the regulation of 
aflatoxin genes. Interestingly, we found one additional 
CRE site in the promoter of mycelial catl in A. parasiticus 
compared with that in the promoter in A. flavus. 

The formation of DNA-protein complexes that 
included AtfB in the aflatoxin and antioxidant promoters 
correlated strongly with the onset of transcription and the 
pattern of transcript accumulation in the target genes. 
These and previous data indicate that AtfB and AflR 
together function as activators of aflatoxin gene transcrip- 
tion (Roze et al. 2004, 2011) and that the interaction of 
AtfB with the target promoters depends on the presence 
of at least one CRE element. 

Our data also strongly suggest that homo and/or hete- 
rodimerization with other transcription factors (AP-1, 
SrrA, AtfA, or other bZIP transcription factors) contrib- 
utes to ability of AtfB to bind the CRE site. Based on the 
data presented in Table 3, ATF-2 (a homolog of AtfA in 
A. nidulans) forms a heterodimer with Jun (a homolog of 
Yapl in S. cerevisiae) or it forms a homodimer at CRE 
sites while ATF-1 (a homolog of AtfB in A. nidulans) 
binds as a homodimer in mammalian cells (Abate et al. 
1990; Hai and Curran 1991; Deppmann et al. 2006). On 
the other hand, Lara-Rojas and coworkers suggested that 
AtfB might interact with AtfA to regulate common target 
genes as AtfB in A. nidulans is similar to Perl in S. pombe 
and Atfl forms a heterodimer with Perl in S. pombe 
(Lara-Rojas et al. 2011). It was shown that Atfl and Perl 
can form a heterodimer and function as a transcription 
activator or repressor for chromatin remodeling in 
S. pombe (Davidson et al. 2004). Thus, these data and our 
previous findings (Roze et al. 2004) strongly suggest that 
AtfB may interact with AtfA in A. parasiticus in response 
to oxidative stress. To address heterodimer/homodimer 
formation between AtfB and AtfA, we will analyze the 
interaction between these two proteins on regulation of 
aflatoxin and antioxidant genes in future work. 

Our EMSA data suggest that Mn sod in A. parasiticus is 
activated by AtfB later than mycelial catl, in agreement 
with findings by others (Kawasaki et al. 1997; Estruch 
2000; O'Brien et al. 2004; Gessler et al. 2007). Mn SOD is 
localized to the mitochondrial matrix and is required for 
detoxification of superoxide generated in mitochondria 



during respiration. We also hypothesize that MsnA may 
assist AtfB binding to the mycelial catl promoter stimu- 
lating its early expression. We wiU analyze this possibility 
in future studies. 

We previously suggested that AP-1 binds at an AP-1 
recognition site 5'-TGAGTAC-3' in the nor-1 promoter 
and together with AtfB activates nor-1 expression (Roze 
et al. 2011). Although this study does not find AP-1 as 
part of a stable DNA-protein complex at the promoters 
of aflatoxin and antioxidant genes, the data do not rule 
out the possibility that AP-1 is a transient component of 
complexes at the target promoters. In yeast, the transcrip- 
tion factor Yapl (Papl ortholog) regulates expression of 
the cytosolic catalase gene CTTl (cttl) via association 
with Skn7 (Prrl ortholog) in response to H2O2 (Lee et al. 
1999; Nguyen et al. 2000; Toone et al. 2001; Aguirre et al. 
2006). Based on this precedent in yeast, A. parasiticus 
AP-1 may function indirectly through a physical associa- 
tion with the ortholog of SrrA {A. nidulans ortholog of 
Skn7 and Prrl in yeasts) to activate mycelial catl and 
aflatoxin biosynthetic genes. In support of this model, 
promoters of the aflatoxin and antioxidant genes (Table 2 
and Fig. 9A) contain AP-1 -motifs located in close prox- 
imity to the SrrA recognition sequence similar to the situ- 
ation in S. cerevisiae (He and Fassler 2005). Alternatively, 
A. parasiticus AP-1 may regulate AflR, a positive regulator 
of the aflatoxin biosynthetic pathway, in order to activate 
expression of aflatoxin genes (Reverberi et al. 2010). We 
will analyze the interaction between AP-1 and SrrA and 
between AP-1 and AflR on the promoters of the aflatoxin 
and antioxidant genes in future studies. 

The temporal formation of DNA-protein complex in 
the Mn sod promoter differed from that observed in three 
other promoters analyzed. By what mechanism does 
A. parasiticus regulate the timing of DNA-protein com- 
plex formation on different target gene promoters? In 
theory, the dynamics of DNA-protein complex formation 
within a specific promoter could be controlled by the 
landscape of specific cfs-acting elements in the promoter. 
In the promoter regions of all five genes analyzed in this 
study, we observed multiple-binding sites for four tran- 
scription factors (AtfB, SrrA, AP-1, and MsnA) located in 
close proximity to each other and within 500 bp of the 
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Figure 9. Proposed model for transcriptional activation of secondary 
metabolism and oxidative stress-response genes by binding of 
transcription factors in response to oxidative stress. (A) Schematic 
representation of the transcription factor-binding sites in the 
promoters of the aflatoxin biosynthetic and antioxidant genes used in 
this study. (B) Regulatory network for transcriptional activation of 
secondary metabolism and oxidative stress-response genes for cellular 
defense against oxidative stress. Based on available experimental 
evidence, we propose that increased levels of intracellular ROS in 
fungal cells downregulate the cAMP-PKA signaling pathway. This 
promotes MsnA binding to stress-response element (STRE) sites in 
promoters of oxidative stress-response (OR) genes including 
antioxidant genes for their activation. Simultaneously, ROS 
upregulates the SAPK signaling cascades. This promotes AtfB and 
SrrA binding (SrrA recruits AP-1) to corresponding CRE, SRRA, and 
API sites in promoters of the oxidative stress-response (OR) genes 
including the antioxidant genes for their induction. Then, MsnA, AtfB, 
and SrrA bind (SrrA recruits AP-1) to corresponding STRE, CRE, SRRA, 
and API sites in promoters of secondary metabolism genes including 
aflatoxin genes for their activation due to excess ROS. Secondary 
metabolites produced in response to ROS play a role in a defense 
mechanism of fungal cells against oxidative stress. PKA, protein kinase 
A; SAPK, stress-activated protein kinase; OR, oxidative stress response; 
SM, secondary metabolism. 



ATG translation initiation codon. Each of the five pro- 
moter regions contained 3-6 AtfB recognition sites, one 
SrrA site, one or two MsnA sites, and multiple AP-f sites. 
The relative position of these specific as-acting elements 
was similar but unique in each promoter (Fig. 9A). In 
three of five promoters regions, CRE sites were detected 
adjacent to ATG (the exceptions were fas-1 and Mn sod); 
in four of five promoters analyzed, an MsnA site was 
detected adjacent (within 100 bp) to ATG (the exception 
was the nor-1 promoter). 

We previously proposed that AtfB aids AflR binding to 
its recognition site in the promoter (Roze et al. 2004) 
and canonical AflR-binding sites 5'-TCGSWNNSCGR-3' 
are present in the nor-1 and ver-1 promoters (Fig. 9A). 
Intriguingly, the A. parasiticus mycelial catl promoter 
contained five motifs that differ from the canonical AflR- 
binding site by only one nucleotide in the 3' end 
(TCGAGTCGCCA, TCGCCAGACGG, TCGTAGCCCCA, 
TCGCAGTGGGA, and TCGTGTTGCCA). The Mn sod 
promoter contained two such motifs (TCGTAGGTGGA 
and TCGTGCGCCAA), and the fas-21 fas-1 intergenic 
region possessed one motif TCGAACACCTA. Currently, 
the functional significance of these motifs is unknown. 
Further analyses are required to determine the binding 
specificity of AtfB to these promoters in vitro and in vivo, 
and also how the landscape of transcription factor recog- 
nition sites affects interaction of AtfB with SrrA, AP-1, 
MsnA, and AflR, its interaction with the promoter, and 
the transcriptional impact during the response to oxida- 
tive stress. 

Aflatoxin biosynthesis is triggered as a part of the cellu- 
lar response to oxidative stress. What is the role of afla- 
toxin biosynthesis in this response? Does it protect the 
cell from oxidative stress? In plants and fungi, secondary 
metabolites play a protective role in adaptation to various 
biotic and abiotic stressors in the environment. In partic- 
ular, plants produce many compounds such as phenyla- 
mides, polyphenol, and flavonoids including anthocyanin 
and quercetin, which have ROS-scavenging properties, for 
their defense (Edreva et al. 2008). Filamentous fungi also 
produce secondary metabolites with antioxidant function 
(Aguirre et al. 2006). For example, the addition of resve- 
ratrol, a natural antioxidant, to Aspergillus ochraceus in 
culture and catalase to Pusarium graminearum in culture 
reduced synthesis of ochratoxin A and deoxynivalenol, 
respectively (Fanelli et al. 2004; Ponts et al. 2007). ROS is 
generated in A. parasiticus under aflatoxin-inducing con- 
ditions early during growth (8-10 h for O^^ generation 
and 18 h for H2O2 generation) (Reverberi et al. 2008). 
Our current data suggest that the expression of mycelial 
catl is activated by MsnA binding at a STRE site and 
AtfB binding at CRE sites in the promoter at 24 h. Based 
on this observation, we propose that antioxidant genes 
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including mycelial catl are activated as the first line of 
defense against ROS formation, and aflatoxin genes are 
activated at a later time point as the second line of 
defense to excessive levels of ROS in agreement with find- 
ings presented by Reverberi et al. (2006). 

Our data obtained in this study and the observations 
of others (Toone and Jones 1998; Thevelein and de Winde 
1999) enable us to expand our previously proposed 
regulatory model (Roze et al. 2011) for regulation of 
stress response and secondary metabolism. We propose 
that increased levels of intracellular ROS activate antioxi- 
dant genes via two signal transduction pathways: First, 
ROS downregulates the cAMP-PKA pathway, which 
results in MsnA binding to the target promoter and acti- 
vation of the antioxidant gene. Simultaneously, ROS 
upregulates SAPK signaling cascade(s), which results in 
AtfB and SrrA binding to the promoters and induction of 
the antioxidant genes. MsnA, AtfB, and SrrA (SrrA 
recruits AP-1) then bind to the promoters of aflatoxin 
biosynthetic genes to assist in their induction (Fig. 9B). 

In conclusion, in fungi, cellular response to oxidative 
stress relies on a complex and dynamic network of tran- 
scription factors that bind in a highly coordinated man- 
ner to the promoters of the target genes and regulate 
their expression. Understanding the mechanisms that 
underlie the temporal regulation of a transcription factor 
binding to its recognition sites in the promoter is one key 
to establishing the link between intracellular ROS forma- 
tion, cellular antioxidant defense, and the benefits for 
survival offered by secondary metabolism. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. Western blot analysis of A. parasiticus protein 
extracts using anti-AP-1 antibodies (peptide 103). Cell 
proteins were extracted from A. parasiticus SU-1 grown 
in YES at 30°C. Western blot analysis was conducted as 
described in Experimental Procedures. Total cell protein 
extracts (90 /ig per lane) and enriched nuclear protein 
extracts (200 /.ig per lane) were resolved by electrophore- 
sis (12% SDS-PAGE). The PVDF membrane with trans- 
ferred proteins was developed with anti-AP-1 antibodies 



(103; 5 fig mL^') as a primary antibody and LI-COR 
goat anti-rabbit secondary antibody conjugated to the 
fluorescent tag IRDye 680 (Li-Cor Biosciences, Lincoln, 
NE). CE, cell protein extracts; EE, enriched nuclear 
protein extracts. 

Figure S2. EMSA of AP-1 binding at the fas-1, ver-1, 
mycelial catl, and Mn sod promoters. Aspergillus parasiti- 
cus SU-1 was grown for 24 or 60 h at 30°C in YES 
medium. Enriched nuclear protein extracts were prepared 
as described in Experimental Procedures. Five micrograms 
of enriched nuclear protein extracts were added to P^^ 
labeled promoter probes of aflatoxin or antioxidant genes. 
Anti-AP-1 antibodies (peptide 103) or preimmune serum 
were added to determine whether these could block 
protein/DNA interaction. (A) Fasl-a probe. (B) Verl-b 
probe. (C) Mcatl-a probe. (D) Msod probe. 
Figure S3. Competition EMSA of fas-1, ver-1, mycelial 
catl, and Mn sod promoters using a 55 bp Verl m pro- 
moter fragment as a competitor. Aspergillus parasiticus 
SU-1 was grown for 24 or 60 h at 30°C in YES medium. 
The 55 bp non-labeled Verl m promoter fragment (con- 
tains STRE site; 50 and 250-fold molar excess) was added 
to compete for labeled Fasl-a, Verl-b, Mcatl-a, or Msod 
probes. (A) Fasl-a probe. (B) Verl-b probe. (C) Mcatl-a 
probe. (D) Msod probe. 

Figure 84. Alignment of mycelial catl and Mn sod pro- 
moter sequences in A. parasiticus and A. flavus. The ATG 
start codon in each sequence is shown with a shadow and 
the CRE sites are underlined. (A) mycelial catl promoter 
sequences. (B) Mn sod promoter sequences. 
Table 81. Primer sequences used to generate double- 
stranded DNA fragment probes for EMSA. 
Table 82. Primer sequences used for qRT-PCR analysis. 
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